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Abstract. The effect of an applied magnetic field on the temperature 
at the maximum of the ZFC magnetization, Mzfc, is studied using the re- 
cently obtained analytic results of Coffey et al. (Phys. Rev. Lett., 80 (1998) 
5655) for the prefactor of the Ncel relaxation time which allow one to pre- 
cisely calculate the prefactor in the Neel-Brown model and thus the blocking 
temperature as a function of the coefficients of the Taylor series expansion 
of the magnetocrystalline anisotropy. The present calculations indicate that 
even a precise determination of the prefactor in the Neel-Brown theory, which 
always predicts a monotonic decrease of the relaxation time with increasing 
field, is insufficient to explain the effect of an applied magnetic field on the 
temperature at the maximum of the ZFC magnetization. On the other hand, 
we find that the non linear field-dependence of the magnetization along with 
the magnetocrystalline anisotropy appears to be of crucial importance to the 
existence of this maximum. PACS numbers : 05.40 +j; 75.50 Tt 



1 Introduction 

Experimental results obtained a few years ago for ferrofiuids |lj and recently for 7-ES2O3 nanoparticles Q 
indicate that for dilute samples (weak interparticle interactions) , the temperature T max at the maximum 
of the zero-field-cooled magnetization Mzfc, first increases with increasing field, attains a maximum and 
then decreases. More recently, additional experiments performed on the 7-E32O3 particles dispersed in a 
polymer ||] confirm the previous result for dilute samples and show that, on the contrary, for concentrated 
samples (strong interparticle interactions) T max is a monotonic decreasing function of the magnetic field 
0j (see Fig. [l]). The behaviour observed for dilute samples (isolated nanoparticles) is rather unusual since 
one intuitively expects the applied field to lower the energy barrier and thus the blocking temperature 
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Fig. 1: The temperature T max at the maximum of the ZFC magnetization plotted against the applied field for 
different sample concentrations. The volume fraction of the particles in the sample 4D (7 — Fe^O^, mean diameter 
~ 8.32 nm) determined from the density measurements is C v = 0.0075, 0.012, 0.047. 



of all particles and thereby the temperature T max Q Resonant tunneling is one of the arguments which 
has been advanced Q] as a mechanism responsible for an increase of the relaxation rate around zero 
field. More recently Mzfc measurements Q have shown an anomaly in the temperature range 40-60K, 
which is probably too high for quantum effects to manifest themselves. Yet another explanation j| of 
the T max effect was proposed using arguments based on the particle anisotropy-barrier distribution. It 
was suggested that for randomly oriented particles of a uniform size, and for small values of the field, the 
field depresses the energy barriers, and thereby enlarges the barrier distribution, so lowering the blocking 
temperature. It was also suggested that the increase of the barrier-distribution width overcompensates 
for the decrease of the energy barrier in a single particle. However, the discussion of the relaxation time 
was based on the original Arrhenius approach of Neel. Here the escape rate, that is the inverse relaxation 
time, is modelled as an attempt frequency multiplied by the Boltzmann factor of the barrier height, with 
the inverse of the attempt frequency being of the order of 1CP 10 s, thus precluding any discussion of the 
effect of the applied field and the damping parameter on the prefactor, and so their influence on the 
blocking temperature. 

It is the purpose of this paper to calculate the blocking temperature using the Kramers theory of the 
escape of particles over potential barriers as adapted to magnetic relaxation by Brown || . The advantage 
of using this theory is that it allows one to precisely calculate the prefactor as a function of the applied 
field and the damping parameter (provided that interactions between particles are neglected). Thus the 
behaviour of the blocking temperature as a function of these parameters may be determined. It appears 
from the results that even such a precise calculation is unable to explain the maximum in the blocking 
temperature T B as a function of the applied field. Thus an explanation of this effect is not possible using 
the Neel-Brown model for a single particle as that model invariably predicts a monotonic decrease of Tg 
as a function of the applied field. 

In view of this null result we demonstrate that the T max effect may be explained by considering an 
assembly of non-interacting particles having a volume distribution. This is accomplished by using Gittle- 
man's Q model which consists of writing the zero-field cooled magnetization of the assembly as a sum of 
two contributions, one from the blocked magnetic moments and the other from the superparamagnetic 
ones, with the crucial assumption that the superparamagnetic contribution is given by a non-linear func- 
tion (Langevin function) of the applied magnetic field and temperature. If this is done even the simple 
Neel-Brown expression for the relaxation time leads to a maximum in T max for a wide range of values of 

* The temperature T max at the maximum of the ZFC magnetization is assumed to be roughly given by the average of 
the blocking temperature Tg of all particles in the (dilute) assembly. 
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the anisotropy constant K. It was claimed in |7| that a simple volume distribution, together with a Neel 
expression for the relaxation time, leads to the same result for FeC particles, although the author used 
the Langevin function for the superparamagnetic contribution to magnetization. 

Therefore, the particular expression for the single-particle relaxation time which is used appears not 
to be of a crucial importance in the context of the calculation of the blocking temperature. 

In the next section we briefly review Kramers' theory of the escape rate. 

2 Kramers' escape rate theory 

The simple Arrhenius calculation of reaction rates for an assembly of mechanical particles undergoing 
translational Brownian motion, in the presence of a potential barrier, was much improved upon by 
Kramers Q . He showed, by using the theory of the Brownian motion, how the prefactor of the reaction 
rate, as a function of the damping parameter and the shape of the potential well, could be calculated 
from the underlying probability-density diffusion equation in phase space, which for Brownian motion 
is the Fokker-Planck equation (FPE). He obtained (by linearizing the FPE about the potential barrier) 
explicit results for the escape rate for intermediate-to-high values of the damping parameter and also 
for very small values of that parameter. Subsequently, a number of authors (0, jl(|) showed how this 
approach could be extended to give formulae for the reaction rate which are valid for all values of the 
damping parameter. These calculations have been reviewed by Hanggi et al.|TT|. 

The above reaction-rate calculations pertain to an assembly of mechanical particles of mass m moving 
along the rc-axis so that the Hamiltonian of a typical particle is 

H=£- + V{q), (2.1) 

where q and p are the position and momentum of a particle; and V{q) is the potential in which the assembly 
of particles resides, the interaction of an individual particle with its surroundings is then modelled by the 
Langevin equation 

* + qP + = A(t) (2 ' 2) 

where X(t) is the Gaussian white noise, and <; is the viscous drag coefficient arising from the surroundings 
of the particle. 

The Kramers theory was first adapte d to the thermal rotational motion of the magnetic moment 
(where the Hamiltonian, unlike that of Eq.(pj|), is effectively the Gibbs free energy) by Brown H in order 
to improve upon Neel's concept of the superparamagnetic relaxation process (which implicitly assumes 
discrete orientations of the magnetic: moment and which does not take into account the gyromagnetic 
nature of the system). Brown in his first explicit calculation || of the escape rate confined himself to 
axially symmetric (functions of the latitude only) potentials of the magneto-crystalline anisotropy so that 
the calculation of the relaxation rate is governed (for potential-barrier height significantly greater than 
ksT) by the smallest non-vanishing eigenvalue of a one-dimensional Fokker-Planck equation. Thus the 
rate obtained is valid for all values of the damping parameter. As a consequence of this very particular 
result, the analogy with the Kramers theory for mechanical particles only becomes fully apparent when 
an attempt is made to treat non axially symmetric potentials of the magneto-crystalline anisotropy which 
are functions of both the latitude and the longitude. In this context, Brown Q succeeded in giving a 
formula for the escape rate for magnetic moments of single-domain particles, in the intermediate-to- high 
(IHD) damping limit, which is the analogue of the Kramers IHD formula for mechanical particles. In 
his second 1979 calculation ^ Brown only considered this case. Later Klik and Gunther [Q, by using 
the theory of first-passage times, obtained a formula for the escape rate which is the analogue of the 
Kramers result for very low damping. All these (asymptotic) formulae which apply to a general non- 
axially symmetric potential, were calculated explicitly for the case of a uniform magnetic field applied at 
an arbitrary angle to the anisotropy axis of the particle by Geoghegan et al. p3[ and compare favorably 
with the exact reaction rate given by the smallest non- vanishing eigenvalue of the FPE |l4| , |l5| , jl(| and 
with experiments on the relaxation time of single-domain particles [Q~7| . 
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In accordance with the objective stated in the introduction, we shall now use these formulae (as 
specialized to a uniform field applied at an arbitrary angle by Geoghegan et al. |H| and Coffey et al. 
pi] , |l5[ ], p6[) for the calculation of the blocking temperature of a single particle. 

A valuable result following from these calculations will be an explicit demonstration of the breakdown 
of the non-axially symmetric asymptotic formulae at very small departures from axial symmetry, mani- 
festing itself in the form of a spurious increase in T max . Here interpolation formulae joining the axially 
symmetric and non-axially symmetric asymptotes (analogous to the one that joins the oscillatory and 
non-oscillatory solutions of the Schrodinger equation in the WKBJ method |l8| ) must be used in order 
to reproduce the behaviour of the exact reaction rate given by the smallest non-vanishing eigenvalue of 
the FPE, which always predicts a monotonic decrease of T max , as has been demonstrated by Garanin et 
al. iflSfl in the case of a transverse field. 



3 Calculation of the blocking temperature 

Following the work of Coffey et al. cited above, the effect of the applied magnetic field on the blocking 
temperature is studied by extracting Tb from the analytic (asymptotic) expressions for the relaxation- 
time (inverse of the Kramers escape rate) plf , which allow one to evaluate the prefactor as a 
function of the applied field and the dimensionless damping parameter r\ r in the Gilbert-Landau-Lifshitz 
(GLL) equation. For single-domain particles the equation of motion of the unit vector describing the 
magnetization inside the particle is regarded as the Langevin equation of the system (detailed in pc} ) . 
Our discussion of the Neel-Brown model as applied to the problem at hand proceeds as follows : 
In an assembly of ferromagnetic particles with uniaxial anisotropy excluding dipole-dipole interac- 
tions, the ratio of the potential energy vU to the thermal energy fc^T of a particle is described by the 
bistable form 

(3U = -a(e ■ n) 2 - £(e • h) (3.3) 

where — t)/(fcgT), v is the volume of the single-domain particle; a — (3K 3> 1 is the anisotropy (reduced) 
barrier height parameter; K is the anisotropy constant; £ = [3M S H is the external field parameter; e, n, 
and h [h = £/2a) are unit vectors in the direction of the magnetization M, the easy axis, and the magnetic 
field H, respectively. 9 and ip denote the angles between n and e and between n and h, respectively. The 
Neel time, which is the time required for the magnetic moment to surmount the potential barrier given 



by (3.3), is asymptotically related to the smallest nonvanishing eigenvalue Ai (the Kramers' escape rate) 
of the Fokker-Planck equation, by means of the expression r w 2tjv / \ 1 where the diffusion time is 



/3M S 



1 

1- Vr 



, ■ (3.4) 

2-f [rir 

7 is the gyromagnetic ratio, M s the intrinsic magnetization, r\ the phenomenological damping constant, 
and r\ r the GLL damping parameter r\ r = 777 M s . 

As indicated above, Brown || at first derived a formula for Ai, for an arbitrary axially symmetric 
bistable potential having minima at 9 — (0,7r) separated by a maximum at 9 rni which when applied to 



Eq.(3.3) for h||n, i.e. a magnetic field parallel to the easy axis, leads to the form given by Aharoni [|2l), 
9 m = cos _1 (-/i), 

Al a ^=a 3 / 2 (I - h 2 ) x [(1 + h) e -«( 1 +' 1 ) 2 + (1 - h) e-^ 1 -^ 2 ] (3.5) 

where < h < 1, h = 1 being the critical value at which the bistable nature of the potential disappears. 

In order to describe the non-axially symmetric asymptotic behaviour, let us denote by (3AU- the 
smaller reduced barrier height of the two constituting escape from the left or the right of a bistable 
potential. Then for very low damping, i.e. for r) r x f3AU- <C I (with of course the reduced barrier 
height /3AU- ^> 1, depending on the size of the nanoparticle studied) we have ||, @, |0| the following 
asymptotic expression for the Neel relaxation time 

rvln - ^- (3-6) 
* Jfix P(Uo - U 1 )e-^-^+^ 2 x P(U - U 2 )e-W°- U ^} 
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For the intermediate-to-high damping, where r/ r x /3AU- > 1 (again with the reduced barrier height 
(3AU- much greater than unity) we have |l6| the asymptotic expression 



where 



IHD 



-P(Ua-Ui) 



-0(Uo-U 2 ) 



(3.7) 
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Here u>\,u)2 and wo are respectively the well and saddle angular frequencies associated with the bistable 
potential, fin is the damped saddle angular frequency and the a? are the coefficients of the truncated (at 
the second order in the direction cosines) Taylor series expansion of the crystalline anisotropy and external 
field potential at the wells of the bistable potential denoted by 1 and 2 and at the saddle point denoted 
by 0. A full discussion of the application of these general formulae to the particular potential, which 

(i) 

involves the numerical solution of a quartic equation in order to determine the c) with the exception of 



the particular field angle ip — j or ^, in Eq^Ey}) is given in Refs. |13j], j22| 

The blocking temperature Tb is defined as the temperature at which r 
time. Therefore, using Eqs.(|3.4| 

T = T n 



T m , T rn being the measuring 



( p . 5| ) and (3.6) (or (3^4), (3J3) and (3/7)) and solving the equation 
for the blocking temperature Tb, we obtain the variation of Tb as a function of the applied field, 
for an arbitrary angle ip between the easy axis and the applied magnetic field. 

In particular, for very small values of ip we have used Eq. ( |3.5|), a s the problem then becomes almost 
axially symmetric and the arguments leading to Eqs.([P|) and (|3.7D fail [||, p3| , ]2(i|] , and appropriate 
connection formulae must be used so that they may attain the axially symmetric limit. We then sum over 
ip, as the easy axes of the particles in the assembly are assumed to be randomly distributed. In Fig. || we 
have plotted the resulting Tb vs. H for different values of the damping parameter r\ r . We have checked 
that lowering (or raising) the value of the measuring time r TO shifts the curve Tb{H) only very slightly 
upwards (or downwards) while leaving the qualitative behaviour unaltered. In order to compare our 
analytical results with those of experiments on particle assemblies, we have calculated the temperature 
T max at the maximum of the ZFC magnetization. In the present calculations we have assumed that M s 
is independent of temperature. We find that the temperature T max behaves in the same way as was 
observed experimentally , || for dilute samples (see Fig. g where the parameters are those of the most 
dilute sample in Fig. [I], with r\ r = 2.5). Moreover, our calculations for a single particle show that the 
blocking temperature Tb (H) exhibits a bell-like shape in the case of intcrmcdiate-to-high damping. This 
behaviour is however spurious as is shown below. 



4 Spurious behaviour of the blocking temperature at low fields 

We have mentioned that initially the non-axially symmetric asymptotic formulae appear to render the 
low-field behaviour of T max . However, this apparent behaviour is spurious as the asymptotic formulae 
(as one may verify (i) by exact calculation of the smallest non- vanishing eigenvalue of the Fokker-Planck 
equation, and (ii) e.g. the IHD formula does not continue to the axially symmetric asymptote) fail at low 
fields, since the IHD formula diverges like h~ 1 / 2 , for all angles, where h is the reduced field as defined in 
sect. |^. The effect of the divergence is thus to produce a spurious maximum in T max as a function of the 
applied field. 

In order to verify this, we have also performed such calculations |ll| , juj , using exact numerical 
diagonalization of the Fokker-Planck matrix. The smallest non-vanishing eigenvalue Ai thus obtained 
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Fig. 2: The blocking temperature Tb as a function of the applied field as extracted from the formulae for the 
relaxation time of a single-domain particle and summed over the arbitrary angle ip, plotted for different values 
of the damping parameter r\ r (see text), and mean volume (V) = 265 nm 3 . (K = 1.25 x 10 5 erg/cm 3 , 7 = 2.10 7 
S^G- 1 , M s = 300 emu/cm 3 , r m = 100 s). 



leads to a blocking temperature which agrees with that rendered by the asymptotic formulae with the 
all important exception of IHD at very low fields where the exact calculation invariably predicts a mono- 



tonic decrease in the blocking temperature rather than the peak predicted by the IHD formula (3.7), so 
indicating that the theoretical peak is an artefact of the asymptotic expansion, caused by using Eq. ( ft.7| ) 
in a region beyond its range of validity, that is in a region where the potential is almost axially symmetric 
due to the smallness of the applied field which gives rise to a spurious discontinuity between the axially 
and non axially symmetric asymptotic formulae. 

An explanation of this behaviour follows (see also |l9|): in the non-axially symmetric IHD asymptote 



Eq.(3.7) which is formulated in terms of the Kramers escape rate, as the field tends to zero, for high 



damping, the saddle angular frequency ujq tends to zero. Thus the saddle becomes infinitely wide and so 



the escape rate predicted by J^q.(3J) diverges leading to an apparent rise in the blocking temperature 
until the field reaches a value sufficiently high to allow the exponential in the Arrhcnius terms to take over. 
When this occurs the blocking temperature decreases again in accordance with the expected behaviour. 
This is the field range where one would expect the non-axially asymptote to work well. 

In reality, as demonstrated by the exact numerical calculations of the smallest non vanishing eigen- 
value of the Fokkcr-Planck matrix, the small field behaviour is not as predicted by the asymptotic 



behaviour of Eq.(3.7) (it is rather given by the axially-symmetric asymptote) because the saddle is lim- 
ited in size to lo . Thus the true escape rate cannot diverge, and the apparent discontinuity between 
the axially-symmetric and non axially-symmetric results is spurious, leading to apparent rise in Tb- In 
reality, the prefactor in Eq. (|3.7| ) can never overcome the exponential decrease embodied in the Arrhenius 
factor. Garanin [p3| (see |19|) has discovered bridging formulae which provide continuity between the 



axially-symmetric Eq.(3.5) and non axially symmetric asymptotes leading to a monotonic decrease of the 
blocking temperature with the field in accordance with the numerical calculations of the lowest eigenvalue 
of the Fokker-Planck equation. 

An illustration of this was given in Ref. Jl9| for the particular case of ip = that is a transverse 
applied field. If the escape rate is written in the form 

t- 1 = —Aexp(-/3AU) 

TT 

where u)% is the attempt frequency and is given by 

2K ^ f\ — w 
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Fig. 3: The temperature T max at the maximum of the ZFC magnetization for an assembly of particles with 
randomly distributed easy axes, as a function of the applied field obtained by averaging the blocking temperature 
in Fig. 2 over the volume log-normal distribution [a — 0.85), see text; with rj r = 2.5. 



then the factor A, as predicted by the IHD formula, behaves as rj r /^/H for h -C l,rjf, while for h = 0, 
A behaves as 2irri r ^/a/i:, which is obviously discontinuous. So that a suitable interpolation formula is 
required. Such a formula (analogous to that used in the WKBJ method |jl8|| ) is obtained by multiplying 
the factor A of the axially symmetric result by e~^Io(^), where 2o(£) is the modified Bessel function of 
the first kind, and £ = 2ah (see Jjl|). 

This interpolation formula, as is obvious from the large and small £ limits, automatically removes the 
undesirable 1 / \[h divergence of the IHD formula and establishes continuity between the axially symmetric 
and non-axially symmetric asymptotes for %j> = 7r/2, as dictated by the exact solution. 

It is apparent from the discussion of this section that the Neel-Brown model for a single particle is 
unable to explain the maximum in T max , as a careful calculation of the asymptotes demonstrates that they 
always predict a monotonic decrease in the blocking temperature. However this effect may be explained by 
considering an assembly of non-interacting particles with a (log-normal) volume distribution and using 
Gittleman's [^) model as shown below, where the superparamagnetic contribution to magnetization is 
taken to be a non-linear function (Langevin function) of the magnetic field. 



5 Possible explanation of the maximum in T max 

Our explanation of the low-field behaviour of T max is based on extracting T max from the zero-field cooled 
magnetization curve assuming a volume distribution of particles. According to Gittleman's model the 
zero-field cooled magnetization of the assembly can be written as a sum of two contributions, one from 
the blocked magnetic moments and the other from the superparamagnetic ones. In addition, we write the 
superparamagnetic contribution as a Langevin function of the applied magnetic field and temperature. 

Gittleman |(| proposed a model in which the alternative susceptibility of an assembly of non- 
interacting particles, with a volume distribution and randomly distributed easy axes, can be written 
as 

oc 

X (T,u) = ±f dVVf(V) Xv (T,Lo), (5.8) 
o 

where Z = J °° dVVf(V), f(V) is the volume distribution function, \v is the susceptibility of the volume 
under consideration, and dVV f (V)xv is the contribution to the total susceptibility corresponding to 
volumes in the range V — V + dV. xv is then calculated by assuming a step function for the magnetic 
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field, i.e. H = for t < and H — Hq — const for t > 0. Then, the contribution to the magnetization 
from particles of volume V is given by 

M v (t) = VH (xo - (xo - XiK t/T ) , (5.9) 

where xo — Mg(T)V/3ksT is the susceptibility at thermodynamic equilibrium and \i = M%(T)V/3Eb 
is the initial susceptibility of particles in the blocked state (see |24j and many references therein). The 
Fourier transform of d5.9| ) leads to the complex susceptibility 

.._fco+fc"Xi) (5J0) 



whose real part reads 



1 + lUJT 



1 + OTT^ 



where r m is the measuring time, and u> is the angular frequency (= 2ttu) 



Starting from (5.11) the application of an alternating field yields: a) x' — Xo if WT "C 1. At high 
temperature the magnetic moments orientate themselves on a great number of occasions during the time 
of a measurement, and thus the susceptibility is the superparamagnetic susceptibility xo- b) x' = Xi 
if lot 3> 1. At low temperature the energy supplied by the field is insufficient to reverse the magnetic 
moments the time of a measurement. Here the susceptibility is the static susceptibility xi - Between these 



two extremes there exists a maximum at the temperature T max . x' can be calculated from (5.11) using the 
formula for the relaxation time r appropriate to the anisotropy symmetry, and considering a particular 
volume V, one can determine the temperature T max . 

In an assembly of particles with a volume distribution, x' can be calculated for a (large) volume 
distribution by postulating that at a given temperature and given measuring time, certain particles are 
in the superparamagnetic state and that the others are in the blocked state. The susceptibility is then 
given by the sum of two contributions 

oo V c 

X '(T, u) = J dVVf(V)xi(T, v) + J dVVf(V) X o(T, u), (5.12) 

V c 

where V c = V C (T, H, v) is the blocking volume defined as the volume for which r = i = r m . x' shows a 
maximum at T max near < Tb > ■ 

If this is done even the simple Neel-Brown[j] expression for the relaxation time leads to a maximum in 
Tmax when the superparamagnetic contribution to magnetization is a Langevin function of the magnetic 
field. Thus the particular expression for the single-particle relaxation time used appears not to be of a 
crucial importance in the context of the calculation of the blocking temperature. 

In Figs. Uj| we plot the result of such calculations (see appendix) where the parameters correspond 
to the samples of Fig. |l]. In Fig. ^ we compare the results from linear and non-linear (Langevin function) 
dependence of the magnetization on the magnetic field. We see that indeed the non-linear dependence on 
H of the superparamagnetic contribution to magnetization leads to a maximum in T max while in the linear 
case the temperature T max is a monotonic function of the field, for all values of K (corresponding to our 
samples) . This only shows that the volume distribution by itself cannot account for the non-monotonic 
behaviour of the temperature T max , contrary to what was claimed in 0. 

In fact, in the non-linear case, T max exhibits three different regimes with field ranges depending on 
all parameters and especially on K (see Figs. ||, ||). For example, for K — 1.5 x 10 5 erg/cm 3 , the field 
ranges are (in Oe) 0.0 < H < 30, 30 < H < 140, and H > 140. In the first range, i.e. very low fields, 
Tmax slightly decreases, then in the second range it starts to increase up to a maximum, and finally for 
very high fields T max decreases down to zero. These three regimes were obtained experimentally in in 
the case of diluted FeC particles. 

^This is the simplest non-trivial case since the relaxation time (and thereby the critical volume) depends on the magnetic 
field. 
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Fig. 4: Temperature T max as a function of the applied field obtained by the calculations of sect. V and app endix . 
Squares : the superparamagnetic magnetization M ap is a lin ear function of the magnetic field given by Eq.(p.2l[). 



Circles : M ap is the Langevin function given by Eq. fl6.22| ). The parameters are the same as in Figs. 1-3, and 
K = 1.0 x 10 5 erg/cm 3 . 



Next, we studied the effect of varying the anisotropy constant K . In Fig. |5| we plot the temperature 
Tmax vs. H for different values of K. It is seen that apart from the obvious shifting of the peak of T max 
to higher fields, this peak broadens as the anisotropy constant K increases. 

We have also varied the volume-distribution width a and the results are shown in Fig. |^. There we 
see that the maximum of T max tends to disappear as a becomes smaller. 

Finally, these results show that the non-monotonic behaviour of T max is mainly due to the non- 
linear dependence of the magnetization as a function of magnetic field, and that the magneto-crystalline 
anisotropy and the volume-distribution width have strong bearing on the variation of the curvature of 
Tmax vs. field. 



6 Conclusion 

Our attempt to explain the experimentally observed maximum in the curve T max (if) for dilute samples 
using the asymptotic formulae for the prefactor of the relaxation rate of a single-domain particle given by 
Coffey et al. [ppf , has led to the conclusion that these asymptotic formulae are not valid for small fields, 
where the maximum occurs. However, this negative result has renewed interest in the long-standing 
problem of finding bridging formulae between non-axially and axially symmetric expressions for the 
prefactor of the escape rate. Recently, this problem has been partially solved in |is[| . 

On the other hand, exact numerical calculations jlj|, H], [^2| of the smallest eigenvalue of the 
Fokker-Planck matrix invariably lead to a monotonic decrease in the blocking temperature (and thereby 
in the temperature T max ) as a function of the magnetic field. We may conclude then that the expression 
of the single-particle relaxation time does not seem to play a crucial role. Indeed, the calculations of sect. 
H have shown that even the simple Neel-Brown expression for the relaxation time leads to a maximum 
in Tmax if one considers an assembly of particles whose magnetization, formulated through Gittleman's 
model, has a superparamagnetic contribution that is a Langevin function of the magnetic field. The 
magneto-crystalline anisotropy and the volume-distribution width have strong influence. 

Another important point, whose study is beyond the scope of this work, is the effect of interparticle 
interactions on the maximum in the temperature T max . As was said in the introduction, this maximum 
disappears in concentrated samples, i.e. in the case of intermediate-to-strong interparticle interactions. A 
recent study based on Monte Carlo simulations of interacting (cobalt) fine particles seems to recover 
this result but does not provide a clear physical interpretation of the effect obtained. In particular, it 
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Fig. 5: Temperature T max (#) for different values of the anisotropy constant K. 
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Fig. 6: T max (_ff) for different values of the volume-distribution width a, and K = 1.5 x 10 5 erg/cm 3 . 

was shown there that interactions have a strong bearing on the effective variation of the average energy 
barrier with field, as represented in an increase of the curvature of the variation of T max with H as the 
packing density (i.e. interparticle interactions) increases. 
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Appendix: Obtaining T max from the ZFC magnetization 

Here we present the (numerical) method of computing the temperature T max at the peak of the zero-field 
cooled magnetization of non-interacting nanoparticles. 
The potential energy for a particle reads 

— = sin 2 9 - 2h cos(?/> - 9) (6.13) 
a 

where all parameters are defined in sect. ||. 

Then, one determines the extrema of the potential U and defines the escape rate A according to 
the symmetry of the problem. Here we consider, for simplicity, the axially-symmetry Neel-Brown model 
where A is given by Eq. (pT 



The next step consists in finding the critical volume V c introduced in Eq. (5.12). V c is defined as the 
volume at which the relaxation time (or the escape rate) is equal to the measuring time r m = 100s (or 
measuring frequency). That is, if one defines the function 

F(V) = \(a,6 a ,9 m ,6 b )-^, (6.14) 

where 9 a ,9 bl 9 m correspond to the two minima and maximum of the potential, respectively, the critical 
volume V c is obtained as the volume that nullifies the function F(V) for given values of T, H and all 
other fixed parameters {^,r] r ,M s and the volume-distribution width a). 
Then, M z f c is defined according to Gittleman's model ||, namely 

V a oo 

Z x M zfc (H, T, ip) = J VV M sp (H, T, V, tl>) + J VV M b (H, T, V, tl>) (6.15) 
v c 

where M sp and Mb are the contributions to magnetization from superparamagnetic particles with volume 
V <V C and particles still in the blocked state with volume V > V c . f(V) = (l/crV2n) exp(— log 2 (V/V rn )/2a 2 ), 
is the lo g-no rmal volume distribution, V m being the mean volume; Z = VV = Vf(V)dV. 
can be rewritten as 

V c Vc oo 

Z x M zfc = J VVM sp - J VVM b + J VVM b . (6.16) 



Now using, 

M b (H, T, V, = sin 2 V, (6.17) 

M b can be taken outside the integral in the last term above. Thus^J, 

v c 

x M b 



Z x M zfc = J VV (M sp - M b ) + Z 
o 

The final expression of M z f c is obtained by averaging over the angle ip ((sin 2 ip} = |), 

7T/2 Vc 

J dtpwa.il> x jVV (M sp - M b ) + (6.18) 



M zfc = - 



o o 



The expression of M sp varies according to the model used. Chantrell et al. |2(| have given an expression 
which is valid for M s HV/k B T < 1, 



MtVH f , , 1 



M sp (H, T, V, i>) = -j^r- (^cos 2 ^ + - 



(6.19) 

*The reason for doing so is to avoid computing the integral which is numerically inconvenient. 



1 - cos 2 V(l - 

^0 
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with 



h = 1 
I a a 



1(a) 



1(a) = 2 dxe 



Note that upon averaging over ip, the expression in (6.19) reduces to 

M 2 VH 
M S JH,T,V) = 3 

3fc B T 

which is just the limit of the Langevin function for M S HV -C fcs? 1 , i.e. 



M sp {H,T,V) = M S C 



k B T 



Therefore, the expression in ( 6.1 8| ) becomes 

M zfc = \ 



VV (M sp - M b ) + 

OA 



(6.20) 



(6.21) 



(6.22) 



(6.23) 



This is valid only in the case of a relaxation time independent of i/j, as in the Neel-Brown model, which is 
applicable to an assembly of uniformly oriented particles. However, if one wanted to use the expressions 
of the relaxation time given by Coffey et al. and others, where r depends on t he an gle ip, as is the case 
in reality, one should not interchange integrations over ip and V, as is done in (6.18), since V c in general 
depends on r and thereby on ip. 

Therefore, the final expression for M z f c that was used in our calculations for determining the tem- 
perature T max is given by eqs. (3.21), (6.22), (6.23). 
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